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The experimental investigation of quantum phases in optical lattice systems provides major chal- 
lenges. Recently, dynamical generation of double occupancy via modulation of the hopping ampli- 
tude t has been used to characterize the strongly correlated phase of fermionic atoms. Here, we 
want to validate this experimental technique with a theoretical study of the driven Hubbard model 
using analytic methods. We find that conclusive evidence for a Mott phase can be inferred from 
such a measurement, provided that sufficiently low temperatures k B T <C t can be reached. 
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The recent progress in cooling and manipulating cold 
fermionic gases in optical lattices allows us to investigate 
phenomena at ever lower temperatures, where intrigu- 
ing many-body effects generate new (quantum) phases 
which are of much interest in condensed matter systems 
as well. While the high tunability of atoms in optical 
lattices is a landmark advantage of these systems, the 
availability of methods allowing for their proper charac- 
terization lags behind the sophistication reached in solid 
state physics. This lack in available experimental probes 
gains importance, now that an increasingly larger class 
of non-trivial quantum phases are in reach. Their defini- 
tive identification requires the translation of the defining 
properties, usually expressed in terms of thermodynamic 
and transport coefficients, to the probes available for cold 
atomic systems. Recently two experiments addressed the 
question of strongly correlated fermions in an optical lat- 
tice |2| . Schneider et al. |2j concentrate on the com- 
pressibility 0], while Jordens et al. [U have made use 
of a new probe, the dynamically generated double occu- 
pancy (DGDO) proposed by Kollath et al. [4J, in order 
to identify a Mott insulator state of strongly repulsive 
fermionic 40 K atoms. This technique, where the signal 
is the change in double occupancy after dynamical mod- 
ulation of the hopping amplitude t, represents a prime 
example of a measurement without analog in classical 
solid state physics |5fl. Using different analytic methods 
in separate regions of the parameter space, cf. Fig. [IJ 
this Letter aims at a theoretical characterization of the 
DGDO scheme and its suitability to describe the transi- 
tion or crossover to the fermionic Mott-insulator. 

Ultracold fermions in an optical lattice are well de- 
scribed by the Hubbard model 



H- 



-t 



(clc ja +h.c.)+U J2 D i = ~tK+UH int , (1) 



where the spin degrees of freedom are realized by two 
hyperfine states. Here, c\ a are fermionic operators cre- 
ating particles with spin a in a Wannier state at site i, 
Di = n^riii measures the double occupancy at site i, 
and the sum runs over nearest neighbors for all N 
lattice sites. Although the phase diagram for the three 



dimensional system is not precisely known yet, the qual- 
itative form was determined (cf. 0, 0] and references 
therein). In particular, at half-filling on a cubic lattice, 
an anti-ferromagnetic phase below T < T N(3cl < 0.3 T F is 
expected; whereas in strongly frustrated systems, typical 
estimates limit the temperature for (paramagnetic) Mott 
physics to a regime below 0.1 T F (T F denotes the Fermi 
temperature) . 

The two defining properties of a fermionic Mott insu- 
lator are the interaction-induced incompressibility and a 
vanishing Drude- weight in the optical conductivity |(|. 
While the first is related to a gap A[i in the charge 
spectrum and can be used as a probe in the cold atom 
setup as well, the measurement of the second one, rely- 
ing on transport properties, is more difficult to perform. 
From this viewpoint, the implementation of the Hubbard 
Hamiltonian with cold atoms not only allows us to study 
this unsolved model in a controlled environment but also 
poses inherently new questions. In comparing theoretical 
and experimental results, it is of interest to know which 
features in the calculation (which is based on a uniform 
system) survive the translation to the (inhomogeneous) 
experimental setup. We thus characterize the DGDO 
technique in a wide parameter regime and conclude that 
frequency-integrated and bond-averaged quantities cap- 
ture the desired information with the least obstructions 
due to confinement-induced changes of line shapes. 
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FIG. 1: (color online) Regimes where the dynamically gener- 
ated double occupancy probe is analyzed: High temperatures 
(I) are discussed in the atomic limit. The regime at zero 
temperature and low densities (II) involves the solution of a 
two-particle problem; the situation at half filling and T — 
(III) is done within a slave-spin mean-field analysis. 
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Our theoretical analysis of the DGDO technique in- 
volves exactly solvable limiting cases of the Hamilto- 
nian (fl]), namely, the atomic limit (t = 0) and the two- 
particle problem, see Fig. [TJ with results relevant in the 
regimes U,T ^S> t, as well as for filling p <C 1 (we set 
k B = 1). On the other hand, in order to discuss the 
Mott physics within the DGDO technique at low tem- 
peratures, cf. Fig. mill) , we introduce and apply a slave- 
spin mean-field approximation which captures the most 
relevant physics on both the t and U scale. Our results 
for the DGDO signal are summarized in Fig. O 

We first derive the expression for the second-order re- 
sponse function describing the DGDO ST>; exciting the 
system during the time A"r mod through a (isotropic) mod- 
ulation H mo4 — 6tcos(uiT)K of the hopping amplitude, 
the task is to calculate the expectation value 

5D(r) = (V(r)|£ A|^(r)) - (Q| £ A|0), (2) 

i i 

of the double-occupancy operator, where |0) denotes 
the unperturbed ground state. We write the perturbed 
wave function ip(r) in the interaction representation, 
|^(t)> = e^(*frr/n)T T e^[-(i/fi)/: oo dT / H mod (T y )]|0) J 
where T T denotes time ordering. Expanding @ to sec- 
ond order in the hopping modulation St and keeping only 
non-oscillatory terms 5D(t) = 7rSt 2 SV(ui)AT IIlod /2h 2 + 
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FIG. 2: (color online) Characteristic DGDO signal for the 
different regimes (I)- (III) of Fig. [T] (a) For the high temper- 
ature region (I), we find an increase in ST> tot for increasing 
U/T and a saturation value 5V tot — » p 2 at U 3> T. (b) 
For the low density phase at T = 0, we find the same be- 
havior with a saturating response at u = U /U c ^> 1 where 
ST> to t p 2 [l — sin(47rp)/47rp]. The response of the system at 
half-filling and T — (III) changes its behavior significantly 
at u = 1. (c) Spectral DGDO signal 8T>(u>) with peaks at 
huo — A and around hu w max(t/ c , {/) identifying transitions 
between the Gutzwiller and upper Hubbard- and those be- 
tween the lower and higher Hubbard bands, respectively, (d) 
Total weights 523^ accumulated in the two peaks in (c). 
Shown in grey is the fluctuation region \u— 1| < u a where our 
mean-field analysis is not valid. 



osc. + 0(St 3 ), we arrive at the spectral transition rate 
(the time-averaged first-order response vanishes) 



J2(n\SD\n)\(n\K\0)\ 2 5(u;-uj n0 ). (3) 



Here, |n) denote the excited states with energy fiw n o = 
huj n — Hluq. Note, that in the current experiment [l[, 
SD(t) is not oc Ar mod , however Q. We define the inte- 
grated quantity 

6V tot = ^f cLu5V(lu), (4) 

normalized to the number Nz/2 of bonds in the system 
(z = coordination number). In the non-interacting limit 
(U = 0), the eigenstates of ^ are those of K and no 
excitations are induced, hence SD tot vanishes and a finite 
interaction U is required to generate a finite result. In 
the following, we evaluate (|4|) in the regimes (I)— (III) . 

In the atomic limit (t = 0, region I), the calculation of 
the change in the double occupancy reduces to a two-site 
problem and ST> tot is given by the square of the probabil- 
ity to find a singly occupied lattice site; its rapid increase 
signals the crossover from weak to strong correlations. In 
the limit U 3> T, ST> tot saturates at p 2 for p < 1 and at 
(2 - p) 2 for p > 1 [cf. Fig. 12 a)]. The atomic limit, valid 
in the regime U,T ^S> t, offers a good starting point for 
analyzing current experiments [H,13|. 

In the low density limit (p < 1, region II), we analyze 
the DGDO by exactly solving the two-particle problem 
with Hubbard interaction [^, [T3|. One then expects a 
signature at finite frequency of order U , as two fermions 
(in a singlet state) on the same site cost an energy U, 
irrespective of the density. We only consider the d = 1- 
dimensional case here, the extension to higher dimen- 
sions is straightforward. The two-particle wave function 
in the singlet channel, ip ecM (r) — (H^i> &2)scat) with r the 
relative coordinate and k\ and ki the momenta of the 
two atoms, can be obtained by solving the Lippmann- 
Schwinger equation (we ignore the center of mass mo- 
tion). Expanding the solution i/) sc * t (r) in Bloch waves 
"0Bioch(?") = (rjfei, feKiock we find the amplitudes 



■ doch (q, k'\q,k) scat = 5 k ,k> + 



(U/N)f(q,k) 



(5) 



£q,k - £q,k> + If] 

with €k, q — 4i cos(g) cos(fc/2) the bare dispersion and 
f(q,k) = [l+U/4:ti\ sin(fc/2) cos(g)|] _1 the scattering am- 
plitude. Here, q = (k± + /c2)/2 and k = k\ — ki are the 
total and relative momenta of the two atoms. In addi- 
tion, a bound state ip bB (r) with amplitudes 

mo o,(q, k'\q) bs = (2U 3 /e*°)V 2 /(e q , k , ~ < s ) (6) 



and energy e^ s — \JU 2 + [4tcos(q)} 2 is found. Rewriting 
K = (H int — H)/t with H iDt = U6 r ,o, we find for the total 
DGDO ((U per singlet pair {k±, k 2 }, 



SD. 
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the above result originates from process where two par- 
ticles in the scattering state \q, fc) Bcat combine to a re- 
pulsively bound pair \q) bs . To obtain the full result at 
small but finite density p, we have to calculate ST> tot = 
53ki,k a n fci,fc2<®t fe oV fe2 > where n kuk2 is the ground state ex- 
pectation value to find a pair {fci, k^) in the singlet chan- 
nel, see Ref. 0. The final result, cf. Fig. |^b), shows 
again an increasing signal for small U /t and a saturation 
at large U/t, with SV tot -> p 2 [l - s±n{4np)/4xp]. While 
the derivation assumes low density, the result up to p = 1 
are shown for comparison with the other approaches Q. 

We now turn to the low-temperature regime of the 
half- filled Hubbard model (JTJ in three dimensions, the 
region III in Fig. [Tj Depending on the amount of Fermi- 
surface nesting, an anti-ferromagnetic phase intervenes 
with the Mott physics. For the optical lattice imple- 
mentation this would require an initial cooling to about 
T{ rs 0.02 T F 0]. Within the slave-spin mean-field treat- 
ment used here, the possibility for the appearance of such 
a magnetic order is ignored. All our results below can 
be obtained within the four-boson approach of Kotliar 
and Ruckenstein [ll| which builds on earlier slave-boson 
methods 

mini- 

However, we apply a minimal formal- 
ism reminiscent of the slave-spin [14j or the slave-rotor 
formulation [l5| which captures all aspects relevant to 
our present discussion. 

On each lattice site, we introduce an auxiliary 
pseudospin-1/2 S with eigenstates (of S z ) |+) (|— )) en- 
coding double and empty (singly) occupied sites with 
excitation energies of order U. In addition, auxiliary 
fermionic operators f a are introduced to describe the low- 
energy quasiparticle degrees of freedom. The physical 
creation (annihilation) operators of the original model 
then are represented 2S x fP. The physical 

states in the enlarged Hilbert-space are |e) = |+, 0), 
\d) = /f/Jl+,0), |T) = /J |- 0), and | j) = /J|-,0), 
where |0) is the vacuum of the /-fermions. The projec- 
tion onto the physical subspace is achieved by imposing 
for each site i the constraint Sf + 1/2 = (J2 a f} a fia ~ I) 2 ' 
as a result, H int involves solely pseudospin operators Sf . 

Within our mean-field solution, we assume prod- 
uct states in pseudo spin and fermion degrees of free- 
dom, thereby relaxing the above constraint. Within 
this approximation, the canonical anti-commutation 
relations are preserved on average, ({c ia , c] CT /}) = 

4(S?Sp({f i(T ,fl,}) = S^Sij, where (...) denotes the 
average over mean-field eigenstates. A a consequence, 
the single-particle spectral weight is preserved as long as 
the spin identity (Sf) 2 = 1/4 is respected. 

We obtain two effective mean-field Hamiltonians: The 
fermion problem assumes the form of a non-interacting 
tight-binding Hamiltonian with a hopping amplitude 
renormalized by a factor gij = A(SfSj), with i,j nearest 
neighbors. On the other hand, the pseudospin problem 



reduces to the transverse Ising model 

H TIM = — JijSfSj + / ] S", (8) 

(hj) « 

with the transverse field h = U/2 and the exchange cou- 
pling J i:j = 4tJ2Mflfja) + c - c )- The transverse Ising 
model is a prime example of a system displaying a quan- 
tum critical point at a critical ratio (2h/J) c = (U/J) c , 
separating a magnetically ordered phase from a quan- 
tum paramagnet. In the following, we restrict our anal- 
ysis to translation- invariant states, for which = J = 
— (16/z) J_ 2 ,de£p CT (e) w 2.67 1, and p a {e) is the non- 
interacting density of states per spin. 

Applying a single-site molecular-field approximation to 
Eq. ||5J), the critical ratio is given by (U/J) c = z and we 
recover the result of the Gutzwiller approximation [l6| 
of the Hubbard model. Within the molecular-field ap- 
proximation, the magnetization of the pseudospin points 
along the z-component of the quantization axis rotated 
around the y-axis due to the action of the magnetic field 
h (or interaction U). The angle of rotation is given by 
tan(tp) = \J\ — u 2 /u for u < 1 and tan(</?) = for u > 1, 
where u = U/U c . The fluctuations around this magneti- 
zation are calculated within a spin- wave analysis 17, 18[ 
and we find the gapped pseudospin-wave spectrum 



fiw(k) 



Uc 

2 



a /1 + 7kM 2 , 
\/u 2 + 7kM, 



for u < 1, 
for u > 1, 



(9) 



with 7k = — (1/^) X)i=i cos(fej) and an excitation gap 
A = hu(0). 

The quantum criticality at u — 1 is reflected in the 
softening of the mode ([9]). For u > 1, the jump A/i 
in the chemical potential from hole to particle doping 
amounts to twice the excitation gap, Ap = 2A; the above 
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FIG. 3: The single-particle spectral density (per spin a) 
A a (uj) calculated within the slave-spin theory for different 
interaction strengths u. In the metallic phase, panels (a) and 
(b), the DGDO probe can induce transitions to the upper 
Hubbard band from either the coherent Gutzwiller band at 
huj — or the lower Hubbard band, see (b). In the insulating 
state, only the latter possibility remains, see (d). 
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pseudospin-mode corresponds to the gapped charge exci- 
tation of the Mott insulator R^, 2(| ■ 

To estimate the validity of the molecular-field plus 
spin-wave calculation, we compare fluctuations with the 
magnitude of the order parameter and define u n through 
the condition l/(z + l)E <0li> W&Sf> + ((5S$) 2 }} = 
(Sq) 2 \ u ii. For the cubic lattice, we obtain u a « 0.9, 
hence, for u « f ± O.f the fluctuation induced correc- 
tions to the molecular-field result are important and the 
validity of the above analysis is limited. 

Below, we will see that the DGDO can be qual- 
itatively understood by the structure of the single- 
particle spectral density, which is defined as A a {uS) — 
-Ek ImG ?( k ' w )/ Ar7r - Here ' G *( k ,u) is the retarded 
single-particle Green's function of the c-fermions which 
is calculated in a straightforward way using the Lehmann 
representation. The result is shown in Fig.[3]in the metal- 
lic as well as in the insulating phase. The gapped mode 
found in the transverse Ising model leads to the incoher- 
ent weight around ± max([/ c , U)/2 in the spectral den- 
sity. In the metallic phase, we find the characteristic 
three peak structure with preformed Hubbard bands cen- 
tered at fko ss ±U c /2 and a coherent Gutzwiller band at 
fko w 0. The Gutzwiller band disappears at u = 1 and 
the Hubbard bands touch at h\o — 0. For u> 1 the Hub- 
bard bands assume a constant width of U c /2 and are sep- 
arated by U. The approximate nature of our treatment 
of the spin problem violates (Sf ) 2 = 1/4 and the spectral 
weight fails to be properly normalized. This failure only 
manifests itself in the fluctuation regime, however. 

The nonlinear DGDO signal is calculated from Eq. (J3j> 
by expressing H int and K in terms of the slave-spin op- 
erators; a careful analysis of the consistency in the order 
of expansion is crucial. For u < 1, we find two distinct 
features in 5T>(ui), cf. Fig. [2j The sharp peak at h\o = A 
with weight 5V< t is due to the excitation of a pseudospin- 
wave mode at k — 0. In the fermionic language, this 
corresponds to the process where the system is excited 
from the Gutzwiller band to the upper Hubbard band, cf. 
Fig.[3jb). The peak measures the quasiparticle weight Z, 
which is absent in the insulating state; hence its vanish- 
ing serves as a clear signal for the metal-insulator transi- 
tion. Without the inclusion of quasiparticle interactions, 
the peak is infinitely sharp. The broad peak around 
fko sa max(f/ c , U) with weight 5T>> ot originates from a 
continuum of pairs of pseudospin waves and corresponds 
to the excitation across the gap between the (preformed) 
upper and lower Hubbard bands, cf. Fig. [31(d) ; this fea- 
ture survives the phase transition and 6T>> t saturates at 
u S> 1. The double-peak structure in Fig. [5] is an arti- 
fact of our mean-field approach [l8} . To obtain the exact 
line-shape a more complicated calculation is required. 

In conclusion, we have calculated the DGDO in various 
regimes and have investigated the suitability of this probe 
to characterize the Mott transition in a fermionic system. 
To make the closest contact to experiment we find that 



it is convenient to study the integrated weight ST> tot . 

For a characterization of the localized phase our find- 
ings show that the weight 5T>> t from frequencies hw w U 
is not an optimal measure for Mott physics: All regimes, 
low density, high temperatures (single-site problem), as 
well as the half-filled case display qualitative similar re- 
sults. Only a quantitative analysis allows to distinguish 
different regimes: the DGDO signal ST>> t in the interac- 
tion dominated regime saturates at a density dependent 
value which is largest at half filling. It turns out, how- 
ever, that one can uniquely identify the Mott transition 
in a strongly frustrated system by making the distinction 
between 8T>< t and 5T>> t . These two weights characterize 
the transitions between the Gutzwiller and upper Hub- 
bard and those between the lower and higher Hubbard 
bands, respectively. The quantity 5T>< t traces the quasi- 
particle weight Z and captures the presence of the co- 
herent Gutzwiller band in the metal; its vanishing then 
serves as a clear signal to identify the Mott transition. 
The main effect of a trapping potential on the disappear- 
ance of this coherent signal is expected to come from a 
residual signal from the low-density regime in the periph- 
ery of the trap and its influence on the close-by insulator. 
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